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The geometric and electronic structures of Najv, Cuat, and Agjv metal clusters are systematically 
studied based on the density functional theory over a wide range of cluster sizes 2 < N < 75. A 
remarkable similarity is observed between the optimized geometric structures of alkali and noble 
metal clusters over all of the calculated cluster sizes A''. The most stable structures are the same for 
the three different metal clusters for approximately half the cluster sizes N considered in this study. 
Even if the most stable structures are different, the same types of structures are obtained when the 
meta-stable structures are also considered. For all of the three different metal clusters, the cluster 
shapes change in the order of linear, planar, opened, and closed structures with increasing N. This 
structural type transition leads to a deviation from the monotonic increase in the volume with N. 
A remarkable similarity is also observed for the N dependence of the cluster energy E{N) for the 
most stable geometric structures. The amplitude of this energy difference is larger in the two noble 
metal clusters than in the alkali metal cluster. This is attributed to the contribution of d electrons 
to the bonds. The magic number is defined in the framework of total energy calculations for the 
first time. In the case of Najv, a semi-quantitative comparison between the experimental abundance 
spectra (Knight et at, Phys. Rev. Lett. 52, 2141 (1984)) and the total energy calculations is 
carried out. The changing aspect of the Kohn-Sham eigenvalues from Af = 2 to A'^ = 75 is presented 
for the three different metal clusters. The feature of the bulk density of states already appears at 
N = 75 for all of three clusters. With increasing N, the HOMO-LUMO gap clearly exhibits an 
odd-even alternation and converges to 0. Although there is a similarity in the N dependence of the 
HOMO-LUMO gap between the three metal clusters, it is much stronger between the two noble 
metal clusters. The growth aspect of the d band below the Fermi level of the noble metal clusters 
with increasing N is presented. A good correspondence is observed in the d characteristic of the 
electronic states between the cluster composed of 75 atoms and the bulk metal. The similarities 
observed in the N dependence of the geometric structures and E{N)s originate from the similarity 
in that of the electronic structures. 



PACS numbers: 

I. INTRODUCTION 



It is beneficial to study the geometric structures and 
various electronic properties of clusters as an interme- 
diate phase of materials between isolated and condensed 
systems; the combination of quantum chemistry and solid 
state physics will result in further developments in the 
field of materials science. In the basis of the density 
functional theory (DFT)ii^, an ab initio many-body the- 
ory for systems in the ground state, a detailed compar- 
ative study of the most and meta-stable geometric and 
electronic structures of three different metal clusters - 
Na^r, Cujv, and Agjv - is systematically carried out over 
a wide range of 2 < iV < 75. The three types of atoms 
composing those clusters have a common feature in that 
the outermost valence electron is one s electron. The 
changing and converging aspects of the geometric and 
electronic structures of alkali and noble metal clusters 
from a diatomic molecule (dimer) to a condensed system 
have been presented explicitly and compared for the first 



time. A simple review of previous studies on clusters is 
presented below. 

In 1984, a study was carried out based on the Woods- 
Saxon type shell model^ to explain the distinctive peak 
observed in the abundance spectra of a Najv cluster at 
special value of iV (iV = 8, 20, 40, 58, 92, ...) evaluated in 
an experiment. Using this shell model, the electronic 
structures of simple metal clusters based on a simple 
spherical effective potential using phenomenological pa- 
rameters. In this model, the Coulomb interactions be- 
tween electrons are completely neglected. In a later 
study, the anormalous TV-dependent stability of a Na^v 
cluster was studied based on various types of shell models 
such as a harmonic oscillator, with the effect of Coulomb 
interactions being neglected^. 

In addition to simple shell models, various types of jel- 
lium models^i^ii were adopted for N&n clusters. In jel- 
lium models for metal clusters, the equilibrium ionic con- 
figuration of a cluster is averaged and replaced by a uni- 
form or nearly uniform potential with a surface bound- 
ary. In most previous studies, the jellium models were 
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incorporated with the DFT and the total energies were 
evaluated. 

In model-based studies of Najv , the iV-dependent clus- 
ter stability was evaluated based on the shell correction 
energy, curvature of the sum of the electronic eigenvalues, 
and total energy^'^^LS.^^lojJ.dLlM^. For the evalua- 
tion of the TV-dependent total energy, a semi-quantitative 
comparison with the peak intensity I{N) in the experi- 
mental abundance spectre^i^ was also carried out. 

However, cluster studies based on these models dif- 
fer from those based on the first principles calculations, 
which directly relate the geometric structures to the elec- 
tronic structures. Even if the experimental magic number 
for a particular value of is obtained accidentally, these 
models are not based on a close relationship between the 
geometric and electronic structurtes of real Najv clusters. 
Furthermore, using these models, it is impossible to study 
the similarities and differences between the alkali metal 
cluster NaAT and the noble metal clusters Cuat and AgN 
because the s electron in the outermost shell as well as d 
electrons in the inner shell contribute to bonding in the 
latter. 

Recenly, the most and meta-stable structures of Najv, 
Cuat, and Ag^r clusters were determined for TV < 22 from 
first principles calculations based on DF T^^i^^i^^i^°i^^ . In 
this study, an extremely systematic comparative study 
based on DFT was carried out for the three different 
metal clusters up to N = 75. Our study successfully 
reproduced the results of stable geometric structures 
obtained in previous studies; furthermore, it success- 
fully discovered the convergence aspects of the electronic 
structures of clusters with those of bulk metals. 

In section II, the computational method and the ap- 
proximation are described. In section III, for the most 
and meta-stable structures are systematically classified 
and the similarities and differences between the three dif- 
ferent metal clusters are described. In section IV, the N 
dependence of the ground state energy of the most stable 
cluster structure is analyzed. In section V, all aspects of 
the electronic structure of the clusters including the odd- 
even alternation of the HOMO-LUMO gap are described. 
In section VI, the conclusion of this study is presented. 

II. COMPUTATIONAL METHODS 

An exact many-body theory for a ground state system 
DFT— - was employed to study the geometric and elec- 
tronic structures of Na^r, Cutv, and AgAr clusters. In this 
theory, the many-body Schrodinger equation is trans- 
formed into a self-consistent Schrodinger - type equa- 
tion for a virtual one - electron in the reference sys- 
tem. This equation is called the Kohn-Sham equation 
and it is derived from the Hohenberg-Kohn variational 
principles, namely, '^'fp^y^ = 0. In the DFT, various 
physical properties of materials such as the exact ground 
state total energy can be evaluated if the exact exchange- 
correlation energy functional is adopted. In principle, the 



exact ionization energy for a system can be evaluated as 
the highest occupied molecular orbital energy by solving 
the Kohn-Sham equation. The other orbital energies are 
considered to approximately describe the quasiparticle 
energy. 

In this study, the generalized gradient approximation 
(GGA) proposed by Perdew and Wang (PW91)^ is em- 
ployed for the exchange-correlation energy functional. 
The spin polarization of the system was also consid- 
ered. In the local density approximation (LDA), which 
was the first approximation developed for the exchange- 
correlation energy functional, the energy functional for 
a real inhomogeneous system is locally approximated by 
an energy functional of a homogeneous electron liquid. 
The main disadvantage of the LDA is the incomplete 
cancellation of the self-interaction fused into the Hartree 
potential by the exchange potential. This leads to an 
incorrect description of the electronic properties when 
electrons are strongly localized in the system. In fact, 
the GGA was developed to overcome this disadvantage 
of the LDA. It improves upon the approximation of sev- 
eral physical properties such as interatomic distance and 
cohesive energy in many systems. 

In order to reduce the number of plane waves used in 
the calculation of the electronic structure, a pseudopo- 
tential method is employed for an approximation of the 
nucleus, inner core electrons, and valence electrons. An 
unempirical pseudopotential called ultrasoft pseudopo- 
tential, which is known to have good transferability^^, 
was employed. In a manner similar to other methods, 
relativistic correction terms such as the mass velocity 
and the Darwin terms are included in the pseudopoten- 
tials. For Na, Cu, and Ag atoms, 3s^, id^'^As^, and 
4ci^°5s"'^electrons were explicitly treated as the valence 
electrons, respectively. 

A local orbital basis set is commonly adopted within a 
real space representation for the electrons of isolated sys- 
tems such as atoms, molecules, and clusters. However, it 
is difficult to converge the total energy effectively because 
some expertise is required to select the basis function for 
each atom. 

Instead of treating the isolated system as is, the to- 
tal energy can be effectively converged by using a plane 
wave basis set within a reciprocal space representation, 
which is frequently used in periodic systems. Here, iso- 
lated systems are approximated by pseudo-crystal sys- 
tems constructed using arranged isolated systems having 
sufficient space between them. 

In this method, the total energy can be effectively 
converged by simply increasing the cut-off energy of a 
plane wave expansion. The total energy of the cluster 
is approximately evaluated from the total energy of the 
pseudo-crystal system per unit cell. Another advantage 
of employing this method is that it is possible to evaluate 
the total energies and electronic structures of the clusters 
and bulk solids using a single method. 

However, the latter method has two disadvantages in 
the study of isolated systems. The first is the incom- 
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pleteness of removing the interaction between isolated 
systems. The second is the inabiUty to separate the total 
energy into the kinetic energy term of electrons, and po- 
tential terms between electron and electron, electron and 
nucleus, and nucleus and nucleus. This is because the 
Ewald method^i is adopted for the evaluation of the to- 
tal energy. In this method, the total energy is effectively 
converged by separating the Coulomb potential into the 
short- and long-range parts, and then treating the former 
in the real space and the latter in the reciprocal space. 

The edge length of each cubic unit cell of the pseudo- 
crystal systems was set to be 20~30 A for Najv clus- 
ters, and 15-^30 A for Cujv and Ag^v clusters. It was 
confirmed that these lengths are sufficiently large to ne- 
glect the electronic wave vector dependence in these elec- 
tronic structures. The cut-off energy for the plane wave 
expansion was set to be 48.7, 233.7, and 180.7 eV for 
NaAT, CuAT, and AgAr, respectively. Brillouin zone in- 
tegration was carried out only for the F-point in the 
pseudo-crystal system. This fc-point selection is known 
to be the most effective in the calculation of an isolated 
system. The residual minimization scheme, direct inver- 
sion in the iterative subspace algorithm was adopted for 
the effective self-consistent calculation of the electronic 
structure a^^i^^ . The convergence criterion of the total 
energy was set to be within 1 x 10~^ eV. 

Geometric structures of the clusters were optimized 
from several hypothetical initial structures using an 
optimization algorithm called the conjugate gradient 
method^i. As initial structures in the geometry optimiza- 
tion procedure for Na^ and Cuat, the optimized cluster 
structures obtained from empirical atomic pair poten- 
tials such as Lennard- Jones, Morse, Sutton-Chen, etci^, 
were employed along with optimized structures based on 
first principles calculations that were obtained in previ- 
ous studie o^^i^^i^" . 

In addition to these structures, structures ex- 
pected to exist from the experimental photoelectron 
spectra32ii2£i2ii2^iS and those expected to exist in the 
neighborhood of the local minimum points on the po- 
tential energy surface from our experience of the opti- 
mization procedures for these systemsi^ were selected. 
The optimized structures of Na^r and Cun obtained in 
this study were employed as the initial structures of Ag^r 
in the geometry optimization procedure, and the inter- 
atomic distances satisfying the nearest neighbor atomic 
distance ratio between the Na, Cu, and Ag bulk crystals 
were scaled. 

The electronic density of states of the Na, Cu, and Ag 
bulk crystals were calculated for comparison with those of 
the clusters. The same cut-off energies and optimization 
method used for the cluster systems were also used for the 
bulk calculations. The lattice constants were optimized 
for the hep and bcc structures for Na, and the fee struc- 
ture for Cu and Ag. Brillouin zone integration was car- 
ried out for the fc-point meshes generated by Monkhorst- 
Pack scheme^. For the bulk bcc, fee, and hep structures, 
the number of meshes were selected to be 8 x 8 x 8, 8 



X 8 X 8, and 8x8x4, respectively. The convergence 
criteria of the total energy for the bulk bcc, fee, and hep 
structure were set to be within 1.0 x 10^**, 1.0 x 10^"*, 
and 5.0 x 10""^ eV/atom, respectively. 

The nearest neighbor diatomic distances, binding en- 
ergies, and bulk moduli for the dimcrs and bulk crystals 
are presented, in the Appendix. 

The Vienna Ab-initio Simulation Package (VASP) was 
employed in this study^. 



III. CLUSTER GEOMETRIES 

A. dependence of most and meta-stable 
structures 

In general, unlike the case of molecules, clusters 
have various energetically nearly degenerated meta- 
stable structures. Therefore, the meta-stable as well 
as most stable structures must be considered for under- 
standing cluster-related phenomena. 

Figures 1, 2, and 3 show the most and meta-stable 
structures of the NaAr, Cun, and Agjy clusters for 3 < 
TV < 12, 13 < iV < 22, and 34 < iV < 75, respectively. 
In these figures, the meta-stable structures are carefully 
selected to show the structural type transition, as de- 
scribed below. Almost all of the most and meta-stable 
structures of the Na^r, Cu^r, and AgAr clusters are sim- 
ilar. Therefore, to save space, the cluster structures are 
represented by the structure of Cuat in the figure. The 
types of these structures are classified by the notations 
I, II, III, IV, and V. This classification obeys the order 
of stability in the structural type of Cuat. The symbols 
L, P, O, and C after I, II, III, IV, and V indicate linear, 
planar, opened, and closed structures, respectively. For 
a detailed analysis of the correlation between the struc- 
tural type and the energetical stability of the clusters, 
these structural types should be further classified accord- 
ing to the symbols. Both opened and closed structures 
are three-dimensional. The coordination number (CN) is 
a necessary concept for the classification of these three- 
dimensional structures. Although the CN can also be de- 
fined for linear and planar structures, it is not required 
for the classification in this study. An opened structure 
is defined as one without atoms whose CNs are greater 
than or equal to 11. Other three-dimensional structures 
are defined as closed structures. The three values listed 
after the symbols L, P, O, and C represent the relative 
energies of the most stable structures of NaA/, Cuat, and 
AgAT, respectively. A value of 0.000 is assigned to the 
most stable structure, and the relative energies of the 
meta-stable structures are expressed in electron volts. 

We consider the case of = 6 as an example. As 
shown in FIG.l, the structural type of the most stable 
structure in Nae is II, while that in Cue and Agg is I. 
The structural types of the 2nd stable structures I in 
Nae are I and in Cue and Agg, II. The structural type 
of the 3rd stable structures in Cug and Agg is III. How- 
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ever, structure III is not found in Nag. A symbol N 
indicates a structural type that has not been found. The 
two different figures of structure II show the same cluster 
structure viewed from different angles. Such different fig- 
ures for the same structure are also shown for the other 
values of N. Structural type IV at = 10 corresponds 
to Plato's polyhedron. The numerous highly symmet- 
ric structures including Plato's polyhedron arc labeled as 
TETRA, OCTA, ICO, CUBO, and DECA, and they re- 
spectively represent a tetrahedron, octahedron, icosahe- 
dron, cuboctahedron, and decahedron, respectively. As 
shown in FIGS. 2 and 3 the three structures composed of 
blue, red, and silver spheres represent the structures of 
NaAf, Cujv, and AgAr, respectively These structures are 
shown to emphasize the relatively large structural differ- 
ences from the relation of similarity. 

The following conclusions are obtained from the most 
and meta-stable structures of the clusters. The optimized 
structures of the N&n, Cun, and Ag^v clusters are al- 
most identical if the most and meta-stable structures are 
simultaneously considered. In particular, the most stable 
structure of the three clusters for N < 22 can be identi- 
fied with high probability. For A^ = 34 and 38, all of the 
most stable structures of the three clusters are different. 
For A^ = 40, the most stable structures of alkali and no- 
ble metal clusters clearly differ. For A^ = 55, all of the 
three different clusters favor structural type I (icosahc- 
dral structure) for their most stable structure. The 2nd 
and 3rd stable structures at A^ = 55 include an atomic 
vacancy in an icosahedral flamed. For all of the three 
different metal clusters, these meta-stable structures are 
more stable than the higher symmetric structures such 
as a decahedron and cuboctahedron. Generally, it is dif- 
ficult to define an atomic vacancy in a cluster because it 
is a system without periodicity in the structure. How- 
ever, an icosahedral structure can be considered to be 
a structure that is cleaved from a quasicrystal. In fact, 
A^ = 55 is the minimum number that shows an atomic 
vacancy stability in the cluster structure. For A^ > 55, 
with increasing A^, structures of the clusters approach the 
periodic bulk crystal structure. Therefore, it is expected 
that the probability of vacancy formation will increase in 
finite cluster systems. 

An overview of the structural type transition of the 
most stable structures in Na^r, Cuat, and Ag^r clusters 
with increasing A^ is described below. The structures of 
NaAT, Cujv, and Ag^r clusters are stable in a linear struc- 
ture at A^ = 2, planar structure at 3 < A^ < 5 (Na), 6 
(Cu, and Ag), opened structure at 6 (Na), 7 (Cu and Ag) 
< A' < 15 (Na and Cu), 16 (Ag), and closed structure at 
A^ > 16 (Na and Cu), 17 (Ag). Systems comprising for or 
more atoms can possibly assume three-dimensional struc- 
tures such as a tetrahedron and trigonal pyramid. How- 
ever, for all of the three different metal clusters, the most 
stable structures are not three-dimensional at 4 < A^ < 5 
(Na), 6 (Cu, and Ag). At A^ = 4, the 3 dimensional struc- 
tures are not stable for these systems. Instead, a rhom- 
bus structure is stable at A^ = 4. A three-dimensional 



structure such as a tetrahedron is realized in other metal 
clusters such as Mg^^. 

The stable cluster structures obtained in this study can 
be roughly understood from knowledge of the electronic 
states using the spherical jellium model. The main dif- 
ference between the first principles model used in this 
study and the spherical jellium model is the manner in 
which the atomic configuration in the cluster is treated. 
In the spherical jellium model, only the outermost va- 
lence electrons in the each atom composing the cluster 
are considered as valence electrons. The other electrons 
and nuclei are unified into an uniform positive charge dis- 
tribution with a spherical surface boundary. Therefore, 
the valence electrons in the spherical jellium model reside 
in the central force field. As a result, each degenerated 
valence electronic state is specified with a monoangular 
momentum: S, P, D, F, in a manner similar to the case 
of electrons in an atom and protons and neutrons in a nu- 
cleus. The stable structures of the metal cluster obtained 
in this study can be roughly considered as those satisfy- 
ing the cluster shapes followed by the valence electronic 
density distribution that originates from the occupied or- 
bitals in the spherical jellium model. 



B. A'^ dependence of the averaged nearest neighbor 
distance and coordination number 

Figure 4 (a) shows the averaged nearest neighbor dis- 
tance (ANND) of the most stable structures of NaAr, 
Cuat, and AgAr clusters for 2 < Af < 75. ANND is de- 
fined as the sum of the nearest neighbor atomic distances 
divided by the number of bonds in a cluster. As shown 
in the figure, the ANND value of the cluster reaches ap- 
proximately 90% of the bulk value at A^ = 20 for all of 
the three different clusters. However, the ANND does 
not converge to the bulk value at A^ = 75, and the values 
arc approximately 98% of the bulk ones. Over the entire 
range of A^, the ratio of ANND among the three different 
clusters agrees well with that of the bulk crystals. 

Figure 4 (b) shows the averaged coordination number 
(ACN) of the most stable structures of NaAr, Cuat, and 
AgTv clusters over the range of 2 < A^ < 75. ACN is 
defined as the sum of the nearest neighbor coordinated 
atomic number for all atoms composing the cluster di- 
vided by A^. The A^ dependence of ACN for the three 
different clusters is similar; it should be noted that the 
ACN values do not reach 9 (75% of the bulk value of 12) 
even at A^ = 75. This is attributable to the existence of 
a surface in the cluster. 



C. Difference between planar structures of sodium, 
copper, and silver clusters 

As described in section HI. A., the most and meta- 
stable structures of the three different metal clusters 
are quite similar. Here, the cluster structures hav- 



5 



ing the same structural type are compared quantita- 
tively using structure III at iV = 7 as an example. 
FIG. 5 shows the details of the structures, represented 
by the Cuj structure as an example. The three val- 
ues represent the angle or interatomic distance of Nay, 
Cut, and Agr, respectively. These values are expressed 
in degrees or angstroms, respectively. The blue and 
green values in parentheses represent the relative an- 
gle ratios: angle (distance)Ara7/angle(distance)cu7 and 
angle(distance)^g7/angle(distance)c«7 , respectively. 

Structure III in Na7 belongs to the C2h point group, 
while in Cu7 and Ag7, it belongs to the D2h point group; 
all of these structures are 3rd stable structures. These 
structures have a Jahn- Teller deformation^ in the equi- 
lateral hexagon in the D^h point group. As expected, the 
relative angle ratio of Ag7 is closer to 1 than that of Na7. 
This is a typical example that indicates the stronger sim- 
ilarity between the two noble metal clusters. However, 
there are several exceptions in clusters having a larger 
value of TV. This may occur due to the increase in the 
number of degrees of freedom in the atomic positions. 
Therefore, it is not possible to present an oversimplified 
picture of the degree of similarity in these cluster struc- 
tures. 



D. TV dependence of the averaged atomic density 

The averaged atomic volume of cluster v{N) is defined 

as 

^(TV) EE ^7r{R)yN. (1) 

Here, (i?) indicates the an averaged distance between 
each atomic coordination Ri and center of mass Rcai 
in a cluster. In other words, (R) is given as 

{R) = ~ ^cm\. (2) 

To compare the TV dependence of v{N) for the differ- 
ent metal clusters, each v(N) is normalized to that of a 
dimer. i.e., v(2). 

Figure 6 (a) shows the TV dependence of a normal- 
ized volume for the most stable structures of NaAr, Cuat, 
and AgAT for 2 < TV < 22. The TV dependence of the 
normalized volume v(N)/v{2) is very similar in the three 
metal clusters. The trends of the change in the TV depen- 
dence at transition sizes TV between different structural 
types arc indicated using the symbols L, P, O, and C. 
v{N)/v{2) decreases significantly at TV from P to O. Fur- 
ther, significant decreases are observed at TV from O to 
C. 

Figure 6 (b) shows the volumes of the most stable 
structures in the range of 15 < TV < 75. Here, the scale 
of v{N)/v{2) is expanded to show the change clearly. 
In general, for TV > 20, v{N)/v{2) increases monotoni- 
cally with TV for all of the three metal clusters. However, 



it should be noted that the values decrease significantly 
with increasing TV between Cu4o (Ag4o) and CU55 (Agss), 
unlike the case of Na4o and Nass. This can be attributed 
to the effect of d electrons in that they may shrink the 
interatomic distances of the quasi-spherical structure-an 
icosahedron which is the most stable structure at TV = 55. 

Figure 6 (c) shows z;(TV)/'i;(2) of the most stable struc- 
ture (icosahedron) of Na^r for TV = 55, 147, and 309 
obtained in our previous study^. Here, v{N)/v{2) in- 
creases monotonically with TV. FIG. 6 (d) shows the 
melting points Tm of the Najv cluster at TV = 55, 147, 
and 309, as observed by Haberland et al^, to discuss the 
relationship with cluster volumes. decreases mono- 
tonically with increasing TV. This TV dependence of Tm is 
in contrast to the trend observed in the case of the cluster 
volume. Aguado et al. showed that the volume v{N) and 
Tm of the most stable structures of the NaAr cluster ex- 
hibited an opposite trend for TV > 55 from the molecular 
dynamics calculations based on DFT-LDA^. This con- 
clusion supports the existence of a relationship between 
the TV dependence of v{N)/v{2) as evaluated by us and 
that of the Tm observed by Habarland et al. Further, the 
more precise DFT-LDA molecular dynamics calculations 
predict that Tm of Na4o is higher than that of Nasg^. Tm 
of Na4o has not yet been observed experimentally. Our 
calculation results of v{N)/v{2) do not contradict with 
the result of the DFT-LDA molecular dynamics study^ 
with regard to the relationship between Na4o, and Na55. 

If the conclusion about the relation between the clus- 
ter z;(TV) or v{N)/v{2) and T„ in NaAr holds for CuAr 
and AgAT for TV < 55, it is expected that T,„ of CU40 
and Ag4o arc lower than that of CU55 and Agss, respec- 
tively. For TV < 22, the volume of the most stable cluster 
changes significantly with the structural type transitions 
(L-^ P-^ C) described in section HI. A. Therefore, 
for TV < 22, the TV dependence of the cluster Tm is ex- 
pected to reflect this change. 



IV. TV DEPENDENCE OF THE CALCULATED 
CLUSTER ENERGY FOR 2 < TV < 75 AND 00 

A. Overall aspects of the cluster binding energy 

As described in section HI. A., the most stable struc- 
tures of the three metal clusters were searched over the 
range 2 < TV < 75 and the ground state cluster ener- 
gies E{N) were evaluated. In general, it was observed 
that E{N)/N of the metal clusters increased with TV and 
approached i?(oo)/c>o, which corresponds to the value of 
the bulk cohesive energy. Interestingly, E{N) exhibited 
higher or lower values at a particular value of TV. To 
investigate the aspects of E{N) in datail, the TV depen- 
dence of the difference between E{N) and the liquid drop 
model (LDM)ia average {E{N)), defined as 

SEiN) = {EiN)) - EiN), (3) 
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is evaluated. Here, {E{N)) is expressed by a linear com- 
bination of three terms as shown below; 

{E{N)) = ayN + a^N^ + a^N^ . (4) 

The first, second, and third terms denote the volume, 
surface, and curvature energy, respectively. The fitting 
parameters a^, a^, and ac were determined as described 
below. Here, the averaged cluster binding energy per 
atom {Eb{N))/N, defined as 

{Ei,iN))/N EE Eil) ~ {E{N))/N (5) 

must be calculated to evaluate SE{N). 

Qy was uniquely determined from the cohesive energy 
of the bulk crystal that corresponds to {Eb{N))/N at 
N ~ oo. Qs and Qc were determined as coefficients of the 
fitting curve {Ei,{N)) /N to the cluster binding energy 
Eb{N)/N over the range of 2 < TV < 75 by applying the 
least squares method and minimizing those values un- 
der the threshold (^h(l)) = 0. FIG. 7 shows the cluster 
binding energy per atom Ei,{N)/N and the LDM average 
{Eb{N))/N as functions of iV^i Eb{N)/N and the av- 
erage {Eb{N)}/N of Naw, Cuat, and AgAr at TV = 2 - 75 
and oo are presented. N^^ = 0,0.237 and 0.794 corre- 
spond to N — oo, 75, and 2, respectively. 

We now consider the difference between Eb{N)/N and 
the LDM average {Eb{N))/N. With an increasing in N 
from 2 to 75 and then cxd, a significant similarity is ob- 
served in the difference among the Najv, Cuat, and Ag^r 
clusters in which each element has one s electron in the 
outermost shell. For the sake of comparison, the exper- 
imental bulk cohesive energy and binding energy of a 
dimer— 1^1^ are also shown in FIG. 7. The values of 
the bulk cohesive energy and binding energy obtained 
through our calculations and experiments are in good 
agreement. However, in an Ag bulk crystal, the calcu- 
lated cohesive energy does not agree well with the exper- 
imental value. This disagreement may be attributable to 
the incompleteness of the description of the Ag {Z ~ 47) 
atom. As reported in previous studies^, the disagree- 
ment originates from the treatment method of the rela- 
tivistic effect through the pseudopotential, the exchange- 
correlation energy functional - GGA, and the basis set for 
the electrons - plane waves. Although an improvement in 
the description is desired, we expect that the qualitative 
feature of the growth behavior of an Ag cluster from an 
atom to bulk shown in this study will not change. 



B. Detailed analysis of the A'^ dependence of cluster 
binding energy 

In section IV. A., a strong similarity is pointed out 
for Najv, Cujv, and AgAr clusters in the N dependence 
of the difference between the cluster binding energy per 
atom-£;b(7V)/A^- and the LDM average-(£'b(A^))/A^. In 
this section, we analyze this similarity in detail. FIG. 7 



shows that the difference for each of the three clusters 
approaches with increasing N . 

However, it is difficult to compare them for large value 
of N . Therefore, the N dependence of the TV multiplied 
values 5E{N) = -E{N) + {E{N)) are noted for the de- 
tailed analysis. 

Figure 8 (a) and (b) show the TV dependence of SE{N) 
for the three metal clusters over the ranges 1 < TV < 22 
and 15 < TV < 75, respectively. In FIG. 8 (a), a strong 
similarity is observed in the TV dependence of 5E{N) for 
1 < TV < 22. These values indicate an odd-even alter- 
nation in TV. In other words, SE{N) increases from an 
odd TV to the next even TV, and then decreases from an 
even TV to the next odd TV. In many cases, SE{N) is pos- 
itive at even TV and negative at odd TV. In FIG. 8 (b), 
a significant similarity is observed in the TV dependence 
of SE{N) for 15 < TV < 75. SE{N) is larger at = 34 
(Na), 55, and 58 and smaller at TV = 40 (Cu and Ag), 
68, 70, 71, and 75 are compared to those at the neigh- 
borhood TV considered in this study. The amplitudes of 
SE{N) for Cun and Ag^r are much larger than that for 
Na^v. Further, the values of Cuat and AgA^ are closer. 
This is attributable to the effect of d electrons in noble 
metal clusters. 



C. Definition of the cluster magic number and 
identification 

In many previous studies of clusters, an N value that 
gives a special cluster stability has been frequently called 
as the magic number. To understand the TV-dependent 
system stability, it is necessary to evaluate the magic 
number based on first principles calculations, because it 
can be a standard magic number. However, the magic 
number has not yet been defined. Therefore, to under- 
stand the TV-dependent system stability, we must first 
define the magic numbers. 

First, the curvature of the cluster energy for TV- 
A2£;(A^), defined as 

A2E{N)=E{N + 1) + E{N-1)-2E{N), (6) 

is discussed. If A2-B(TV) exhibits a positive peak at TV, TV 
may be a magic number because E(N) might be a local 
minimum for TV. FIG. 8 (c) shows the TV dependence of 
A2E{N) for the three metal clusters. Here, the magic 
numbers of the three clusters can be identified as TV = 
2, 4, 6 (Na, Ag), 8, 10, 12, 14, 18, and 20. For these values 
of TV, an odd-even alternation is clearly observed. All 
of these values of TV are even numbers and it is natural 
to consider that the odd-even alternation originates from 
the shell closing of each electronic state in the clusters. 
In addition, the TV dependences of the peak intensities 
are similar in the three clusters, although the absolute 
values of the peaks are larger in the noble metal clusters 
than in the alkali metal cluster for most values of TV. As 
in the case of 6E{N) described in IV. B., the values of 
A2E{N) are closer in the noble metal clusters. 
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Although it is possible to derive other conclusions from 
the absolute values of A2E{N), only A2E{N) may not 
be a sufficient criterion for the evaluation of magic num- 
bers. In general, additional information is required for 
the identification. 

SE{N) can be a criterion for identifying magic numbers 
because of its definition: the energy difference between a 
real cluster and the continuously averaged energy model, 
LDM. The magic numbers N should satisfy the condition 
dE{N) > 0. This condition is satisfied for = 2, 4 (Ag), 
6, 7, 8, 10, 12 (Cu), 14, 17 (Cu), 18, 19, 20, and 21 (Na). 

If only A2E{N) and 6E{N) are considered, magic 
numbers N can be defined as 



A2E{N) > 0,SE{N) > 0. 



(7) 



If this definition is employed for the magic numbers in 
Na^v, Cutv, and Ag^r clusters for N < 21, the common 
magic numbers A^ = 2, 4 (Ag), 6, 8, 10, 12, 14, 18, and 20 
are identified. 

In this study, A2E{N) for A^ > 34 is not evaluated. For 
N > 34, magic numbers are identified only from SE(N). 
For these values of A^, we consider those that satisfy 



SE{N) > 



(8) 



to be magic numbers. Therefore, A^ = 34 (Na), 38 (Na), 
40 (Na), 55, and 58 are identified as magic numbers. In 
this definition, the odd number A^ = 55 is also considered 
to be a magic number. 

In our first principles calculations, the common feature 
of the magic numbers in Najv, Cu^v, and Agjv can be at- 
tributed to the delocalized s valence electrons. The local- 
ized d electrons in the noble metal clusters increase the 
stability difference between magic and not-magic clus- 
ters. 



D. Af-dependent stability of Najv clusters from the 
experiment and total energy calculations 

The A^ dependence of the peak intensity I{N) in 
the experimentally observed abundance spectra, and the 
cluster energy E{N) evaluated from the first princi- 
ples calculations have not yet been compared quantita- 
tively. This comparison is important for understanding 
the magic number observed in the experiment. On the 
basis of several assumptions, semi-quantitative compar- 
ison methodi^ii^ has already been developed for an ex- 
perimental result of Najv"^ and for the theoretical total 
energy calculations. Although experimental results for 



Cu^v and Ag 



44.45 



are available, they cannot be eas- 



ily compared with our calculation results quantitatively. 
Therefore, in this section, we only focus on the magic 
numbers of Na^f. 

With regard to experiments with Najv, only the result 
of Knight et al^ is discussed. As shown in other exper- 
iments with Najv such as those by Bj0rnholm et alJ^, 
Rabinovitch et aZ4^, etc., main peaks were also observed 



at A' = 8,20,40,58, and 92. Therefore, without signif- 
icant improvements in the experimental method, these 
abundance spectra and magic numbers appear to remain 
unchanged. Obviously, if the experimental results remain 
unchanged, the theory itself must be reconsidered. 

Here, the curvature of E{N) evaluated from our cal- 
culations and the peak intensity I{N) in the abundance 
spectra were scmi-quantitatively compared based on the 
method described below. The abundance spectra of Na^v 
clusters reported by Knight et al. are considered to be 
attributed as follows. First, Na atomic vapor is formed 
by heating Na bulk solid in an Ar-gas-fiUed closed space 
in the experimental apparatus^. Then, NaA? clusters 
arc formed by the adiabatic cooling of the mixture after 
passing through the skimmer in the apparatus. Then, A^- 
distributcd neutral charged Najv clusters were assumed 
to be obtained in a thermal equilibrium condition. If the 
assumption for the NaAr distribution is true, the relation- 
ship 



p{N) = Aexp 



keT ) 



(9) 



must be satisfied. Here, the p{N), A, ks, and T denote 
the number density of neutral charged Na^v, a constant, 
Boltzmann constant, and the absolute temperature, re- 
spectively. These clusters were ionized by light and then 
accelerated by an electric field for mass selection. In these 
steps, the A^-distribution of p{N) is assumed to remain 
unchanged. Finally, these clusters were detected and the 
A^-distribution of the neutral charged clusters was ob- 
served using the detector. Therefore, finally. 



liN) = Bp{N) 



(10) 



holds. Here, i? is a constant. From equations (6), (9), 
and (10), the following equation is derived. 



A2E{N) = keTln 



I{NY 



I{N+1)I{N - 1) 



(11) 



Namely, A2E{N) can be evaluated from I(N) of the ex- 
perimental abundance spectra. 

Figures 8 (d) and (e) show the A^ dependence of 
A2E{N) evaluated experimentally by Knight et al- 
A2EExpt.{N) and by the first principles calculations in 
this sindy-A2EFP-DFT{N) for A" < 22 and 15 < A^ < 
75, respectively. For the latter discussion, 5E{N) evalu- 
ated from the first principles calculations in this study, 
SEfp^dft{N), and from the DFT-based spherical jel- 
lium model (SJM) calculations by Genzken et alJ^, 
5Esjm-dft{,N), are also shown in these figures. 

We first compare the A^ dependence of A2EExpt.{-^) 
and A2Efp-dft{N). Generally, there is a good agree- 
ment between them. However, contrasting trends are 
observed in two rows of A^: (1) A^ = 5,6 and 7 and (2) 
A^ = 17 and 18. The peak intensity of A2EExpt.{-^) ^^'^ 
A2Efp-dft{N) is relatively strong at A^ = 8 and 20. 
For A^ > 34, A2Efp-dft{N) is not evaluated in this 
study. 



8 



As shown in Figures 8 (d) and (e), there is a similar- 
ity in the N dependence between /S.2Epp_jjpt{N) and 
5Epp-.dft{N). Considering this similarity, instead of 
^2Epp-dft{N), 5Epp-p,ft{N) can be considered as 
a comparative value to l^2EExpt.{N). 

In a manner similar to 1^2Epp~dft{N)^ 
5Efp-dft{N) exhibits good agreement with 
/S.2EExpt.{-^) with regard to the N dependence. How- 
ever, in a manner similar to the case of A2Efp-dft{N), 
6Efp-dft{N) exhibits a different N dependence as 
compared to /S.2EExpt.{^) ^'^o rows of N: (1) 
iV = 5 and 6 and (2) = 17 and 18. Here, it should 
be noted that the order of peak intensity at A = 8 
and 20 of 5Efp~dft{N) exhibits better agreement 
with A2l?_Ea;pt.(A) than with /S.2Efp-dft{N). Based 
on these results, 5Efp~dft{N) is compared with 
/^2EExptXN) for A > 34. 

In /S.2EExpt.{N)^ distinctive peaks arc formed at 
A = 40 and 58. Similar to the case of A2EExpt.{N), 
6Efp-dft{N) is relatively large at A = 58. However, 
5Efp-dft{N) at a = 40 is not particularly so distinc- 
tive. Although A2i?_Expi.(A) is remarkably distinctive 
at A = 40, the value of 5Efp~dft{^^) is similar to 
that of 5EFP-DFT{i^)- Further, for A = 34 and 55, 
5Efp-dft{N) exhibits distinctive peaks and these val- 
ues are larger than SEfp-dft{4:0) and SEfp-dft{58)- 
In addition, for A = 68,70,71, and 75, 5Efp-dft{N) 
has negative values, and these features are not exhibited 
by A2EExpt.{N). 

As shown here for Najv, although the A dependence 
of 5Efp-dft{N) exhibits good agreement with that of 
A.2EExpt.{-^) in the range of A < 22, the same is not 
necessarily true in the range of A > 34. This point is 
discussed in the next subsection. 



E. Improvement of the magic number description 
for Najv clusters from first principles calculations 

Wc evaluate the source of this disagreement in the A 
dependence of 5E{N){/S.2E{N)) obtained from theoreti- 
cal total energy calculations and A2E{N) obtained from 
the experiment. In this subsection, we review the A- 
dependent system stability of NaAr by the improvement 
of theoretical models. 

As shown in Figures 8 (d) and (e), 5Esjm-dft{N) 
obtained by Genzken et alJ^ exhibits an A dependence 
similar to that of SEfp-^dft{N) for 2 < A < 58. How- 
ever, the A dependence oi 5Efp-dft{N) exhibits a bet- 
ter agreement with that of A2EExpt.{N) for A = 34,40, 
and 58 as compared to that of 5Esjm-dft{N). This 
resut suggests that the consideration of the explicit ionic 
configuration in the total energy calculation is significant 
to realize a better agreement with the experiment with 
regard to the A-dependent system stability of Najv- 

Figures 8 (f) and (g) show the A dependence of 
A2E{N) for NaAT evaluated from various DFT-bascd 
models such as the spherical jellium model (SJM-DFT) 



by Chou et al^, spheroidal jellium model (S'JM-DFT) 
by Ekardt et al^, uhimate jellium model (UJM-DFT) 
by Koskinen et al^, and the first principles model (FP- 
DFT) in this study for 1 < A < 22 and 15 < A < 
75, respectively. These A2E{N) values are denoted as 

A2EsjM-DFt{N), A2Es',JM-DFt{N), A2EUJM-DFT, 

and A2Efp~dft{N), respectively. The various jel- 
lium models are differ from the first principles calcula- 
tions in the treatment of the ionic configuration and the 
exchange-correlation energy. 

As shown in Figures 8 (f) and (g), all of the calcula- 
tion results exhibit distinctive peaks at A = 8, 18, and 
20. In these figures, remarkable differences are obtained 
in the number of peaks between A2E{N) of the spheri- 
cal models, SJM-DFT, and non-spherical models, S'JM- 
DFT, UJM-DFT, and FP-DFT. From these figures, it 
is apparent that the improvement in the A dependence 
in A2E{N) from the jellium models to the first princi- 
ples model that explicitly treats the ionic configuration 
is smaller as compared to that from the spherical models 
to the non-spherical models. 

In a manner similar to the case of 5Efp-dft{N) and 
A2Efp~dft{N) for 1 < A < 22, as shown in Figures 
8 (c) and (g), the A dependence of A2Esjm-dft{N) 
is similar to that of 5Esjm-dft{N) over the range of 
2 < A < 58. Here, it is assumed that there is a similarity 
between the A dependence of 5Efp-dft{N) and that 
of A2Efp-dft{N) for 34 < A < 58. 

As shown in Figure 8 (e), A2Esjm-dft{^0) 
{A2Es',jm-dft{^0)) is positive. However, it is 
smaller than A2Esjm^dft{'20), A2Esjm-dft{34), 
and A2£'s,/a/-_d_ft(58) (A2£'s' ja/-_d_ft(20), 
A2Es'jm-dft{M), and A2ES' jm-dft{5S)). On 
the other hand, as shown in FIG. 8(g), 6Esjm~dft{N) 
for A = 40 is negative. Therefore, based on the defini- 
tion of the magic number given in equation (7), A = 40 
is not a magic number in SJM-DFT. However, A = 40 
may be a magic number in S'JM-DFT and FP-DFT 
although SEfp-dft{4:0) has a small positive value. 
Although the magic feature at A = 40 is expected to 
be weak in S'JM-DFT and FP-DFT, it is also expected 
that the consideration of the deviation from central 
force field in the total energy calculation is significant 
for obtaining a better agreement in A-dependent system 
stability with the experiment. 

Here, it should be noted that A2i?s,/M-nFT(68) is 
nearly eV. As shown in FIG. 8 (e), there is a 
large difference between A2EExpt.{^^) eV and 
6Efp-dft{Q8) ^ —0.5 eV. From the similarity be- 
tween the A dependence of SE{N) and A2E{N), 
A2-E_fp-_d_ft(68) is expected to differ significantly from 
A2-E_Ea;pt. (68) as compared to A2£'sjM-_D_Fr (68). There- 
fore, it can be said that the A dependences of 
A2Efp-dft{N) and 5Efp~dft{N) are not necessar- 
ily more similar to that of A2EExpt. (A) as compared to 
those of A2Esjm-dft{N) and 6Esjm-dft{N). 

We now consider the possibilities for reducing the 
difference between the theoretically and experimentally 



9 



obtained TV-dcpcndent system stability. First, we dis- 
cuss the possibility of improving has the magic num- 
ber description based on the evaluation of the ground 
state energy. In the case of Nrat, it has already been 
shown that the N dependence of the experimental abun- 
dance peak at A'^ = 5,6, and 7 cannot be expressed 
even in terms of A2E{N) of the configuration interac- 
tion calculations^Si^. Therefore, it is expected that the 
improvement of the ground state energy evaluation will 
never lead to and improvement in the agreement for small 
values of N. 

Thus for, several researchers have shown stable struc- 
tures of Na4(j^i^i^. In previous studies, a nearly spher- 
ical structure with high symmetry has not yet been ob- 
tained as the most stable structure. Instead of such a 
structure, we have found a structure with low symme- 
try (III) that is 0.320 eV more stable than structure (II), 
which is similar to the most stable structure found thus 
far—. Since the energy differences between these struc- 
tures arc small, the most stable structure. III, did not 
exhibit a strong peak for iV = 40 in 5Efp-dft{N). 
Therefore, we expect that 5Epp-oft{^^^) will not real- 
ize a strong magic feature even if a more stable structure 
is found in Na4o. 

Within SJM, the N dependence of the energy per va- 
lence electron for Na^r was evaluated from a more pre- 
cise first principles calculation method called the diffu- 
sion Monte Carlo (DMC) calculation^. A weak magic 
feature was obtained for = 40 as a dip in the total 
energy vs. N curve. From the result a weak magic fea- 
ture is expected for = 40 in the N dependences of 
A2EsjM-DMc{N) and 5Esjm-dmc{N), in a manner 
similar to the cases of SJM-DFT and FP-DFT. How- 
ever, by considering the explicit geometrical structure in 
DMC, the magic feature for = 40 will be strengthened 
as compared to that for A^ = 34 and 58. 

Second, we discuss the possibility of realizing improve- 
ments by a more proper evaluation of the system stability. 
In the experiments, the temperature in the apparatus is 
expected to be related to the final results. However, in 
the cluster study based on the total energy calculations, 
the temperature of the system is K. If the considera- 
tion of the temperature is critical for the A^-dependent 
system stability, instead of the internal energy E{N), the 
free energy F{N) must be considered. Further, the most 
stable and meta-stable structures of Na^r must both be 
considered properly for the A^-depcndent system stabil- 
ity. To consider the contribution of mcta-stablc struc- 
tures, information about the potential energy surface or 
free energy landscape is required. 

Finally, we discuss the possibility of realizing improve- 
ments from other viewpoints. In the case of Najv, various 
shell models such as the Woods-Saxon type^iiii and har- 
monic oscillator type^ that neglect the Coulomb interac- 
tions between electrons exhibit a magic feature at A^ = 40 
in the A^ dependence of the shell correction energy and 
the curvature for the sum of the electronic eigenvalues. 
These A^ dependences are strongly affected by the A^ de- 



pendence of the energy gap between the highest occupied 
molecular orbital (HOMO) and the lowest unoccupied 
molecular orbital (LUMO), i.e., the HOMO-LUMO gap. 
As shown in FIG. 9, within FP-DFT, the HOMO-LUMO 
gap for the most stable structure of Na4o is larger than 
that of Na34. The absolute value of the HOMO-LUMO 
gap may be larger if the quasi-particlc energy evaluated 
from GWA is considered^. The improvement is signifi- 
cant if I{N) in the abundance spectra of Najy is a value 
that is more strongly related to the shell correction en- 
ergy and the curvature of the sum of electronic eigenval- 
ues as compared to A2E{N) or SE{N) evaluated from 
E{N). 

As described in this section, by a more proper treat- 
ment of the system within the total energy calculations, 
discrepancies in the description of the A^-dependent sys- 
tem stability with the experimental results in the descrip- 
tion of can be reduced. However, this is difficult even in 
the case of the simplest metal cluster Najv shown here. 
For discussing the A^-dependent cluster stability, various 
problems must first be solved. However, the magic num- 
ber for the ground state, which avoids empirical param- 
eter, is significant because it can be a standard to un- 
derstand the A^-dependent system stability. Further, for 
practical reasons, the evaluation method based on first 
principles calculations is expected to retain the value in 
predictions of the experimental values of magic numbers. 



V. ELECTRONIC STRUCTURES 

A. Evolution of the electronic structure from atom 
to bulk 

As described before, remarkable similarities are ob- 
served in the A^ dependence of the most and meta-stable 
structures and in the ground state energies of Najv, Cuat, 
and AgAT clusters. These similarities arc attributed to 
the A^ dependence of the electronic structure. In this 
section, A'^-dependent electronic structures for the most 
stable structures of each metal cluster are discussed. 

Figures 9 (a), (b), and (c) show the Kohn-Sham en- 
ergy spectra for the most stable structures of the three 
metal clusters over the range of 1 < A^ < 75 and oo 
(bulk). The density of states of the bulk crystal obtained 
from the band calculation is shown in the right-hand side 
space of each figure for the sake of comparison. Each 
Fermi level of the bulk crystal is set to the next HOMO- 
LUMO gap for the N = 75 cluster. For all of the three 
figures, the space between HOMO and LUMO is colored 
blue. The red lines represent the occupied and unoccu- 
pied electronic energy levels. 

For most values of A^, the HOMO-LUMO gaps are 
small at odd values of N and large at even values of A^, 
which results in the odd-even alternation of the HOMO- 
LUMO gaps. However, the HOMO-LUMO gap is large 
for several odd values of A^ and small for even values 
of A^. It should be noted that the odd-even alternation 
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of the HOMO-LUMO gaps in Cun and Ag^r are sig- 
nificantly larger than that of Na^r for the benefit of d 
electrons. Further, the similarity in the odd-even alter- 
nation is more distinctive between the noble metal clus- 
ters. Generally, the HOMO-LUMO gap of each metal 
cluster converges to the Fermi level of the bulk crystal 
with increasing N. 

For all of the metal clusters, with increasing N, the 
feature of each bulk energy band gradually appears in the 
electronic states. The energy width between the bottom 
and top of the occupied energy levels almost converges to 
that of the bulk crystal at iV = 75. In the case of Najv, 
over a wide energy range, the states arc characterized by 
s and p. In Cujv and Agjv, in addition to the s, and 
p characterized states similar to those shown in Na^v, 
energetically localized d band type states also appeared. 
It should be noted that the features of localized d band 
in the bulk crystals already appeared at a rather small 
value of N in the width of the d states and the energetical 
distance from the top of the d states to the Fermi level. 
With increasing N, the strongly d characterized states 
gradually expand and almost converge to the width of 
the d band in each bulk crystal at iV = 75. 

For the same type of geometrical structures as those 
described in section II. A., distinctive simiralities were 
also observed in the underwent splitting manner in which 
electronic states that are strongly characterized by s and 
p. As shown in FIG. 9, the energetical distances for 
the adjacent states strongly characterized by s and p are 
larger in NaAr than those of Cu^v and Ag^r. This rela- 
tionship can be roughly understood from the relationship 
between the energetical distance for the adjacent states 
and the width of a well in the quantum well model, as 
described in the previous paper—. In this case, the s and 
p characterized electrons correspond to the quantum in 
the well. As shown in FIGS. 4 (a), 5, and 6 (a) and (b), 
for the existence of spatially and enegetically localized 
d electrons, closer interatomic distances are obtained in 
noble metal clusters as compared to those in alkali metal 
clusters. Namely, the width of the well is narrower in 
the noble metal clusters. As a result, larger energetical 
distances are realized in the noble metal clusters. 

The stability of the spin polarized state of the clusters 
is summarized as given below. For any odd TV clusters, 
the stability of the spin-polarized state is higher than 
that of the spin-unpolarized state. On the other hand, for 
most even N clusters, the energy of the spin-unpolarized 
state is lower than that of the spin-polarized state. How- 
ever, several structures that have high dcnegeracy in the 
neighborhood of the HOMO in the electronic state are ex- 
ceptions. For a geometric structure with high symmetry, 
such as an icosahedron at = 13 and TV = 55, cubocta- 
hedron at TV = 38, capped icosahedron at TV = 71, and 
Marks decahedron at TV = 75, the energy levels around 
HOMO are highly degenerated. Therefore, the high-spin 
states exhibit a higher stability than the low-spin states 
in systems such as an icosahedron for Nasg^ and CU55, 
although the case is opposite for Agss. However, for all of 



them, the energy differences between the high- and low- 
spin states are very small-0.033, 0.022, and 0.012 eV, 
respectively. The system is stabilized by lowering the 
symmetry from 1^ to C'i in the structures based on the 
Jahn- Teller theorem^. 



B. HOMO-LUMO gap, E{N), and v{N)/v{2) for the 
most and meta-stable clusters 

In section V. A., only the electronic structure of the 
most stable cluster structure is discussed. In this sec- 
tion, the electronic structures of the most and meta- 
stable cluster structures are compared to understand the 
relationship between the geometric and electronic struc- 
tures of the clusters. To save space, only the result of 
NaAT is presented. FIG. 10 shows the HOMO-LUMO 
gap ((a) and (b)), A2£;(TV) (a'), .^^^(TV) ((c) and (d)), 
and v{N)/v{2) ((e), and (f)) of the most and meta- 
stable structures of NaAr clusters for 2 < TV < 22 and 
15 < TV < 75, respectively. The values of the most sta- 
ble structures are connected by a line. Several important 
features are observed from the comparison of the most 
and meta-stable structures, as described below. 

A close correlation is observed between the HOMO- 
LUMO gap and A2i?(TV) of the most stable structure. 
As shown in FIGS. 10 (a) and (a'), there is a strong 
correlation in the TV dependence between them for a 
range of 2 < A^ < 21. Although the HOMO-LUMO 
gap considered in this study is the Kohn-Sham HOMO- 
LUMO gap, originally, the HOMO-LUMO gap is a phys- 
ical value that corresponds to the difference between the 
ionization potential and the electron affinity in the sys- 
tem. The ionization potential of a neutral charged clus- 
ter composed of TV atoms is defined as the energy dif- 
ference between the total energy of the neutral cluster 
energy- i? (TV, n)- and the one electron detached charged 
cluster-i?(TV, n — 1). The electron affinity of a neutral 
charged cluster composed of TV atoms is defined as the 
energy difference between the total energy of the neutral 
cluster- £^ (TV, n)- and the one electron attached charged 
cluster-£;(TV, n 1). Therefore, the HOMO-LUMO gap 
of a cluster composed of TV atoms- HLG{N) is defined 
as HLG{N) = E{N,n + 1) + E{N,n - 1) - 2E{N,n). 
On the other hand, the energy curvature of a cluster 
composed of TV atoms-A2i?(TV)- is defined as A2-E(TV) = 
E{N + 1, n) + E{N-l, n)-2E{N, n). Although they are 
clearly different physical values, the forms of these two 
types of values are very similar. The close correlation 
described above should be studied in detail in the future. 

Generally, as shown in Figures 10 (a) and (b), the most 
stable structures exhibit a relatively large HOMO-LUMO 
gap. Further, as shown in FIG. 9, a cluster that has a 
large HOMO-LUMO gap tends to have a lower HOMO 
and higher LUMO in the electronic structures, as de- 
scribed in section V. A. However, from a comparison with 
the energy differences between the most and meta-stable 
structures from SE{N) shown in FIGS. 10 (c) and (d). 
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it is apparent that a large HOMO-LUMO gap is not a 
necessary condition for the most stable structure. For 
example, in Na4, the meta-stable structure III exhibits a 
larger HOMO-LUMO gap than that of the most stable 
structure I. Although it is not shown in this paper, the 
same type of example can also be shown for Cujv and 

AgAT. 

Finally, the N dependence of the normalized cluster 
volume v{N)/v{2) of the most and meta-stable struc- 
tures, and the relation to the total energy and the 
HOMO-LUMO gap is noted. As shown in FIGS. 10 (c) 
and (f), it is apparent that the cluster volume is directly 
dependent on the structural type, namely, L, P, O, and 
C. From the system stability, as shown in FIGS. 10 (c) 
and (d), it is apparent that the most stable structures 
do not necessarily have a minimum value of v{N)/v{2). 
Further, the N dependence of v{N)/v{2) ioi N >7 docs 
not change significantly if energetically closed structural 
isomers to the most stable structures are also considered. 
As shown in the figure, it is apparent that the relation- 
ship between the HOMO-LUMO gap and v{N)/v{2) for 
the most and meta-stable structures cannot be simplified. 

VI. CONCLUSIONS 

The A^-dependent geometric structure, system stabil- 
ity, and electronic structures of NaAr, Cuat, and AgAr 
metal clusters are studied in detail for a range of 2 < A^ < 
75 based on the density functional theory. Strong sim- 
ilarities are observed between the three different metal 
clusters. These similarities originate from the outermost 
is s electron in the alkali and noble metal atoms com- 
posing the each cluster. Much stronger similarities are 
observed between the two noble metal clusters for the 
benefit of d electrons. The most stable structures are 
the same for the three different metal clusters for ap- 
proximately half the cluster sizes A'^ considered in this 
study. Even if the most stable structures arc different, 
the same type of structures are obtained if the meta- 
stable structures are also considered. For all of the three 
clusters, the structural type of the most stable structure 
changes in the order L ^ P — > O ^ C with increas- 
ing A^. This structural type transition leads to a devia- 
tion from the monotonic increase in volume with A^. A 
remarkable similarity is also observed for the A^ depen- 
dence of cluster energy E{N) for the most stable geomet- 
ric structures. This similarity is related to the similarity 
in the electronic structures. The amplitude of this en- 
ergy difference is larger in the two noble metal clusters 
than in the alkali metal cluster. This is attributed to 



the contribution of the d electrons to the bonds. The 
magic number is defined in the framework of total en- 
ergy calclations for the first time. In the case of NaAr, a 
semi-quantitative comparisson between the experimental 
abundance spectrai^. and the total energy calculations is 
carried out. For the improvement of the agreement with 
the experimental result, several possiblities arise for the 
total energy calculations. The changing aspects of the 
Kohn-Sham eigenvalues from A^ = 2 to A^ = 75 are pre- 
sented for the three different metal clusters. The features 
of the bulk density of states already appeared at A^ = 75 
for all of the three clusters. With increasing values of 
N, the HOMO-LUMO gap clearly exhibits an odd-even 
alternation and converges to 0. This alternation is a spe- 
cific feature of alkali and noble metal clusters in which 
each element has one s electron in the outermost shell. It 
is also found that in some cases, the HOMO-LUMO gap 
is large at odd value of A^ and small at even value of A^. 
Although there is similarity in the A^ dependence of the 
HOMO-LUMO gap between the three metal clusters, it 
is much stronger between the two noble metal clsuters. 
The spatially localized d electrons in the noble metals 
are energetically localized. The growth aspect of the d 
band below the Fermi level of the noble metal clusters 
with increasing A^ is presented. A good correspondence 
is obtained in the d characteristic of the electronic states 
between the cluster composed of 75 atoms and the bulk 
metal. 
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and bulk crystals for Na, Cu, and Ag in the equilibrium position evaluated from the DFT— '^GGA(PW91— ) calculations in 
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FIG. 1: The most and meta-stable structures of Najv, Cujv, and Agjv clusters for 3 < A'^ < 12 are shown. Almost all of the 
most and meta-stable structures of these clusters are similar. Therefore, the cluster structures are represented by the structures 
of Cuiv to save space in this figure. The structures are classified according to their types using the notations I, II, III, and IV. 
This classification obeys the stability order in the structural type of Cujv. The symbols L, P, O, and C used after the notations 
indicate linear, planar, opened, and closed structures, respectively. Opened and closed structures are both three-dimensional 
structures. An opened structure is defined as one without any atoms whose coordination number is greater than or equal to 
11. Other three-dimensional structures are difiend as closed structures. The three values following the symbols L, P, O, and 
C represent the relative total energies of the most stable structures of Najv, Cujv, and AgAr, respectively. A value of 0.000 
is assigned to the most stable structure and the energy of meta-stable structures is expressed in units of electron volts. The 
symbol N is used to denote a structure that cannot be identified. The two different figures shown above show views of the same 
cluster from different angles. The numerous highly symmetric structures, including Plato's polyhedron, are labeled as TETRA, 
OCTA, ICO, CUBO, and DECA, and they respectively represents a tetrahedron, octahedron, icosahedron, cuboctahedron, and 
decahedron. 
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FIG. 2: The most and meta-stable structures of Naiv, Cujv, and Agiv clusters for 13 < A'' < 22 are shown. The notations used 
are the same as those in FIG[T] Some isomers of NaAr and Agjv are also shown for cases where their structural deviations are 
relatively large from the relation of similarity. The structures of Najv and Agjv are indicated in blue and silver, respectively. 
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FIG. 3: The most and meta-stable structures of Nslm, Cujv, and AgM clusters for 34 < Af < 75 are shown. The notations used 
are the same as those in FIGS [l] and [2] 
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FIG. 4: (a) shows the A'^ dependence of the averaged nearest neighbor distance (ANND) in units of angstroms, and (b) shows 
the averaged coordination number (ACN) from atom to bulk for Najv, Cujv, and Agjv for 1 < A'^ < 75, and oo (bulk). The 
ACN values of the bulk crystal - 12 (Na: hep, Cu: fee, Ag: fee) and 8 (Na: bee), are not shown in (b). References a, b, and c, 
correspond to ref.—, ref.—, and ref.—, respectively. 
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FIG. 5: The deviation from the relation of similarity shown in structure III of Nay, Cuy, and Ag7 shown in Figure 1. 
The values of the angles and interatomic distances are shown in the order of Nay, Cut, and Agy. Structure III of Cut 
is shown as an example. The angles and interatomic distances are expressed in units of degrees and angstroms, respec- 
tively. The blue and green values in the parentheses represent the relative ratio of angle(distance)jva7/angle(distance)cii7 and 
angle(distance) Agy /angle(distance)cii7 , respectively. 
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FIG. 6: (a) and (b) show the A'^ dependence of the normahzed cluster volume v{N)/v{2) for the most stable structures of Na]v, 
Cujv, and Agjv for 2 < N < 22 and 15 < A'^ < 75, respectively. Each structure is denoted using the symbols L, P, O, and C. 
(c) shows the A'^ dependence of v{N)/v{2) for Na]v (A^ = 55, 147, and 309), while (d) shows the TV dependence of the melting 
point Tm observed by Haberland et al^ 




FIG. 7: The dependence of the binding energy of the ground state structures of Najv, Cujv, and Agjv (1 < Af < 75, oo (bulk)) 
on A'" 3. For each element, the liquid drop model (LDM) average is also shown. References a and b correspond to ref.— and 
ref.— , respectively. 
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FIG. 8: (a) and (b) show the dependence of SE{N) for the most stable structures of Najv, Cujv, and Agjv clusters at 
1 < < 22 and 15 < A*' < 75, respectively, (c) shows A2E{N) for the most stable structures of Najv, Cujv, and Agjv at 

1 < AT < 21. (d) and (e) shows the A^ dependence of A2E{N) = fci3rinjp^:iiy^L_ of Najv as given by Knight et al. (ref. a^) 

at 700 kPa Ar in the cluster production step, 5E{N) of the spherical jellium model (ref. b^^) and first principles model, and 
A2-E(A'') of the first principles model for Najv clusters at 1 < A'' < 22 and 15 < AT < 75, respectively. For A2E{N) value given 
by Knight et al., the value of the temperature is set to 800 K as measured mearured in the nozzle channel before the cooling 
step in the experiment, (f) and (g) show the N dependence of A2i?(A'') of Najv clusters evaluated by various jellium models 
and the first principles model for 1 < A'' < 22 and 15 < A'^ < 75, respectively. References c, d, and e, correspond to ref.—, ref.—, 
and ref.—, respectively. 
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FIG. 9: A'^ dependence of the Kohn-Sham eigenvalues from atom to bulk calculated for the most stable structures of Najv, 
Cujv, and Agjv at 1 < TV < 75 and oo (bulk). The red lines indicate occupied and unoccupied electronic energy levels. For 
each spin-polarized system, the up spin states (left) and down spin states (right) are separated by a thin black line. The space 
between HOMO and LUMO is indicated in blue. Each projected value of s, p, and d to the Kohn-Sham state is connected by 
water, green, and blue lines, respectively. The density of states of each bulk crystal is shown to the right of each figure. Here, 
the Fermi level is represented by a dashed line. The Fermi level of a bulk crystal is arranged near the HOMO of the clusters 
composed of 75 atoms. Each colored line in the bulk crystal is the same as that in the case of clusters. 
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FIG. 10: (a) and (b) show the dependence of the HOMO-LUMO energy gaps of Najv for the most and meta-stable structures 
at 1 < A'^ < 22 and 15 < < 75, respectively, (a') shows the A'' dependence of A2E(N) for the most stable structures of Najv 
at 1 < Af < 21. (c) and (d) show the A'^ dependence of SE{N) for Najv for the most and meta-stable structures at 1 < A'^ < 22 
and 15 < Af < 75, respectively, (e) and (f) show the A'^ dependence of v{N)/v{2) for Najv for the most and meta-stable 
structures at 1 < Af < 22 and 15 < A'^ < 75, respectively. Each different type of structure, namely, linear (L), planar (P), 
opened (O), and closed (C), is indicated in purple, black, red, and blue, respectively. 



